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ABSTRACT 

Magnetic reconnection is commonly considered as a mechanism of solar (eruptive) flares. A deeper 
study of this scenario reveals, however, a number of open issues. Among them is the fundamental 
question, how the magnetic energy is transferred from large, accumulation scales to plasma scales 
where its actual dissipation takes place. In order to investigate this transfer over a broad range 
of scales we address this question by means of a high-resolution MHD simulation. The simulation 
results indicate that the magnetic-energy transfer to small scales is realized via a cascade of consecutive 
smaller and smaller flux-ropes (plasmoids), in analogy with the vortex-tube cascade in (incompressible) 
fluid dynamics. Both tearing and (driven) "fragmenting coalescence" processes are equally important 
for the consecutive fragmentation of the magnetic field (and associated current density) to smaller 
elements. At the later stages a dynamic balance between tearing and coalescence processes reveals a 
steady (power- law) scaling typical for cascading processes. It is shown that cascading reconnection 
also addresses other open issues in solar flare research such as the duality between the regular large- 
scale picture of (eruptive) flares and the observed signatures of fragmented (chaotic) energy release, 
as well as the huge number of accelerated particles. Indeed, spontaneous current-layer fragmentation 
and formation of multiple channelised dissipative/acceleration regions embedded in the current layer 
appears to be intrinsic to the cascading process. The multiple small-scale current sheets may also 
facilitate the acceleration of a large number of particles. The structure, distribution and dynamics of 
the embedded potential acceleration regions in a current layer fragmented by cascading reconnection 
are studied and discussed. 

Subject headings: Sun: flares — Acceleration of particles — Magnetic reconnection — Magnetohydro- 
dynamics (MHD) — Turbulence 



1. INTRODUCTION 

It is generally conjectured that solar flares represent 
a dissipative part of the release of the magnetic en- 
ergy accumulated in active regions at the Sun. The 
'standard' CSHKP model fsee. e.g. IShibata fc Tanumal 
12001b iMandrinil f2010t iMagara et atlll996L and references 
therein) agrees well with the observed large-scale dynam- 
ics of eruptive events. In this model the flares are initi- 
ated by eruption of a flux-rope (in many cases observed 
as a filament ) via, for example, a kink or torus insta- 
bility fe.g. iKliem fc Torokl 12005 |Torok fc Klieml 120051: 
IWilliams et al.l 120051 : iKliem et al.l 12010ft evolving later 
into a Coronal Mass Ejection (CME). The latter is 
trailed by a large-sc ale current layer behind the eject a 
([Lin fc Forbesl l2000ft . In this trailing current layer re- 
connection is supposed to give rise to various observed 
phenomena like hot SXR/EUV flare loops rooted in Ha 
chromospheric ribbons, HXR sources in loop-top and 
foot-points, a nd radio emiss i ons of various t ypes. Some 
authors (e.g. iKo et al.l [2001 iLin et al,ll2QQ7h argue that 
the bright thin ray-like structure observed sometimes be- 
hind CMEs may represent a manifestation of the density 
increase connected with this current sheet (CS). 

However, closer analysis of the classical CSHKP model 
revealed some of its open issues. Namely, the time-scales 
of reconnection in such a thick flare CS appeared to be 
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much longer than typical flare duration. In other words, 
reconnection rate in such configurations has been found 
to be insufficient for the rapid energy release observed 
in flares. Later it was found that the dissipation neces- 
sary for reconnection in the practically collisionless solar 
corona is an es sentially plasma-kinetic process (see, e.g., 
Buchner 2006) which takes place at very small spatial 
scales. Hence, the question arises, how sufficiently thin 
CSs can build up within the global-scale, thick CME- 
tr ailing current layer: Open is the actual physical mech- 
anism that provides the energy transfer from the global 
scales, at which the energy is accumulated to the much 
smaller scales, at which the plasma-kinetic dissipation 
takes place. 

Addressing these questions IShibata fc Tanumal (2001) 
suggested a concept of cascading (or fractal, as they call 
it) reconnection. According to their scenario a cascade of 
non-linear tearing instabilities occurs in the continuously 
stretched current layer formed behind a CME. Multiple 
magnetic islands (helical flux-ropes in 3D), also called 
plasmoids, are formed, interleaved by thin CSs. Due to 
increasing separation of the plasmoids in the continu- 
ously vertically extending trailing part of the CME the 
interleaving CSs are subjected to further filament ation 
until the threshold for secondary tearing instability is 
reached. This process continues further, third and higher 
levels of tearing instabilities take place, until the width 
of the CSs reaches the kinetic scale. 
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This scenario has recently been suppor ted by the ana- 
lytical theory of plasmoid instability by Loureiro et al. 
(2007). They show that the high-Lundquist-number 
systems with high enough current-sheet length-to- 
width ratio are not subjected to the slow Sweet- 
Parker reconnection but they are inherently unstable 
to formation of plasmoids on very short time-s cales. 
Samtanev et al.l (I2009D. iBhattachariee et al.l (|2009l ). and 
Huang fc Bhattacharjed (|2010f ) confirmed predictions of 
this analytical theory by nu merical si mulat ions with 
high Lundquist numbers. iNi et al.l (|2010l ) general- 
izes the model b y presence of shear flo ws around cur- 
rent sheet (CS). llJzdenskv et al.l (|2010l ) relate the the- 
ory of plasmoid instability further to the concept of 
fractal r econnection suggested by IShibata fc Tanumal 
(200lh. iShepherd fc Cassakl (|201Q[ ) and iHuang et all 
(|201Q[ ) study the plasmoid instability numerically at 
smaller scales and investigate its relation to the Hall re- 
connection. They found various regimes of parameters 
where different type of reconnection prevails. 

Eventually, however, kinetic scales are reached 
where dissipation and particle acceleration take place 
most likely via kinetic coalescence of micro-pla s moids 
and, possibly, their shrinkage (iDrake et al.l 120051 : 
iKarlickv fc Bartall2007l : Ikarlickv et al.ll2010D . 

In addition to the issue of energy transport there are 
also other questions that remain open in the CSHKP 
model. It is its apparent insufficiency to accelerate 
such a number of particles in its single diffusion re- 
gion around the X-line that would correspond to the 
fluxes i nferred from HXR observations i n the thick-target 
model (jFletcherl 120051 : iKrucker et all I2008L and refer- 
ences therein). And, furthe rmore, the HXR and radio 
(e.g. decimetric spikes , see IKarlickv et al.l 1 19961 120001 : 
iBarta fc^ Karlickv l200ll ) observations indicate that the 
particle acceleration takes place via multiple concurrent 
small-scale events distributed chaotically in the flare vol- 
ume rather than by a single compact acceleration process 
hosted by a single diffusion region. Such observations are 
usually referred to as signatures of fragmented/chaotic 
energy release in flares. 

Because of these difficulties an alternative concept 
based on the so-calle d "self-organized cr it icality" (SOC) 
has been proposed (|Aschwandenl 120021 : IVlahosI 120071 ). 
This class of models is based on the idea of multiple 
small-scale CSs embedded in chaotic (braided) magnetic 
fields that are formed as a consequence of random mo- 
tions at the system boundary (photosphere). Multiple 
CSs can host multiple reconnection sites what provides 
natural explanation of observed signatures of fragmented 
energy release. At the same moment they provide larger 
total volume of diffusion regions, perhaps sufficient to 
account for the observed particle fluxes also quantita- 
tively. Organised large-scale picture (i.e. coherent struc- 
tures like the flare-loop arcades) should be in the case of 
SOC-based models achieved by the so-called avalanche 
principle: A small-scale energy release event can trigger 
similar events in its vicinity provided the system is in 
marginally stable state (due to the continuous pumping 
of energy and entropy through the boundary). Never- 
theless, it is difficult to achieve such coherent large-scale 
structures as they are usually observed in solar flares in 
the frame of SOC-based models. 

Thus, solar flares appear to be enigmatic phenomena 
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Fig. 1. — Large-scale magnetic reconnection from the point of 
view of theory of dynamical systems. Schematic view of the cascade 
of energy transfer from large to small scales and the window of 
scales resolved in our simulations. Lower abscissa shows the wave 
number k, corresponding characteristic scales A = 2ir/k are at the 
upper abscissa. The numeric values at the scale axis correspond to 
scaling based on typical coronal parameters (for details on scaling 
see Section [2]). 



exhibiting duality between the regular, well-organized 
dynamics of flares observed at large scales and signatures 
of fragmented/chaotic energy release seen in observations 
related to flare-accelerated particles. While the coherent 
global eruption (flare) picture seems to be in agreement 
with the CSHKP scenario, the observed fragmented- 
energy-release signatures favor the SOC-based class of 
models. 

In the present paper, we suggest that cascading re- 
connection can address these three pressing questions 
(i.e. energy transport across the scales, accelerated- 
particles fluxes, and the organized/chaotic picture du- 
ality) as closely related to each other. In our view the 
energy is transferred from large to small scales by the cas- 
cade of fragmentation of originally large-scale magnetic 
structures to smaller elements. We identify two elemen- 
tary processes of this fragmentation (see Section [3]). In 
the course of this process also the initial current layer 
fragments into multiple small-scale, short-living current 
sheets. These current sheets are hierarchically embed- 
ded inside the thick current layer in qualitatively self- 
similar manner. In this sense the cascading fragmenta- 
tion reminds SOC models, but now the chaotic distribu- 
tion of small-scale currents results from internal instabil- 
ities of the global current layer. The fragmented current 
layer represents the modification of the standard CSHKP 
model and thus it keeps coherent large-scale picture of 
solar flares. At the same time it addresses the observed 
signatures of fragmented energy release and the question 
of efficient particle acceleration. We believe that the cas- 
cading reconnection in solar coronal current layers can 
thus address the three main problems mentioned above 
en bloc, and it reconciles the two concepts of the stan- 
dard CSHKP and SOC-based models seen hitherto as 
antagonistic. 

The paper is organized as follows: First, we describe 
the model used in our investigations. Then we present re- 
sults of our high-resolution MHD simulation of cascading 
reconnection in an extended, global, eruption-generated 
current layer. We identify the processes that lead to 
the fragmentation of magnetic and current structures to 
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smaller elements. Then we analyze the resulting scaling 
law of the energy cascade. We describe the structure, 
distribution and dynamics of small dissipation regions 
embedded in an initially thick current layer. Finally, we 
discuss the implications that cascading reconnection have 
for theory of solar flares. 

2. MODEL 

Generally speaking, the solar flare involves three kinds 
of processes that take place in different scale domains 
- see Fig HJ At the largest scales magnetic-field energy 
is accumulated. During this stage flux-rope (filament) 
is formed and its magnetic energy increases. Eventually 
it looses its stability and gets ejected. This process al- 
ready represents (ideal) release of the magnetic-energy 
at large scales. Subsequently, a current layer is formed 
and stretched behind ejected flux-rope. Energy transfer 
from large scales at which the magnetic energy has been 
accumulated to the small dissipation scale occupies in- 
termediate range of scales. The dissipation itself takes 
place at smallest, kinetic scales. 

In this paper we aim at studying energy transfer from 
large to small scales by means of numerical simulations. 
Despite the high spatial resolution our simulation is still 
within the MHD regime. Also, we do not address the 
very process of energy accumulation - i.e. the flux-rope 
formation and energization, nor its instability and sub- 
sequent current-layer formation. Instead we assume a 
relatively thick and extended current layer to be already 
formed at the initial state of our study. In order to cover 
a large range of scales we limit ourselves to the 2D ge- 
ometry allowing for all three components of velocity and 
magnetic field (commonly referred to as 2.5D models). 
This is a reasonable assumption since observations show 
that the typical length of flare arcades along the polarity 
inversion line (PIL) is much larger than the dimension 
across the PIL. 

In the range of scales that we are interested in the 
evolution of magnetized plasma can be adequately de- 
scribed by a se t of compre ssible resistive one-fluid MHD 
equations (e.g. Priest 1984): 
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The set of equations (pQ) is solved by means of Finite 
Volume Method (FVM). For the numerical solution it 
is first rewritten in its conservative form. The (local) 
state of magneto-fluid is then represented by the vector 
of basic variables \l/ = (p, pu, B, J7), where p, u, B, and 
U are the plasma density, plasma velocity, magnetic field 
strength and the total energy density, respectively. The 
energy flux S and auxiliary variables - plasma pressure 
p and current density j - are defined by the formulae: 
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Fig. 2. — Regions of adapt ively enhanced resolution (thick-black- 
line bounded areas) on the background of current density mag- 
nitude (gray-scale) and magnetic field lines (dashed). The right 
panel shows zoomed view of the selected rectangle. Only the rel- 
evant sub-set of whole computational domain is shown, note the 
strongly anisotropic axes-scaling selected in order to show the high- 
resolution sub-domains better. 
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and g is the gravity acceleration at the photospheric 
level. Microphysical (kinetic) effects enter into the large- 
scale dynamics by means of transport coefficients - here 
via a (generalized) resistivity rj. In general, the role 
of non-ideal terms in the generalized Ohm's law in- 
creases as the current density becomes more concen- 
trated via current sheet filamentation. To quantify this 
intensification we use the current-carrier drift velocity 
= \j\/( en e) as the threshold for non-ideal effects to 
take place. Such behavior is presumed by theoretical con- 
siderations and confirmed by kinetic (Vlasov and PIC 
codes] numerical experiment s ([Buchner fc Elkinal 12005. 
120061 : iKarlickv fc Bart a 2008). In particular, we assume 
the followin g; law for (generalized) resistivity (see also 
Klie m et al.feoOOf ): 
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In order to study the energy-transfer cascade it is ap- 
propriate to cover a large range of scales. For struc- 
tured grids it means the utilization of very fine meshes. 
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Fig. 3. — Projection of the initial state (t = 0) to the xz-plane. 
Black lines represent the magnetic field lines, red areas the magni- 
tude of the current density, its scale is given underneath. Enhanced 
view on selected area (the right panel) can be directly compared 
with the third panel in Fig. [4] 



For a given simulation box size the number of finite grid 
cells is limited technically by CPU-time and memory de- 
mands. Alternatively, one can use a refined mesh only 
at locations where the small-scale dynamics becomes im- 
portant - this idea forms the base for th e Adaptive Mesh 
Refin ement (AMR) technique (see, e.g.jBerger fc Oligerl 
1984; iDreher fc Grauerl 12005b Ivan der Hoist fc Keppensl 
2007). We implemented the numerical solver for the 



MHD system of equations (pQ) in the form of block AMR 
Finite- Volume Method (FVM) code: Whenever the cur- 
rent sheet width drops below a certain threshold, a re- 
fined mesh sub-domain is created and initialized with 
values one step backward in time. Its evolution is then 
computed using accordingly refined time-step (this pro- 
cedure is commonly known as sub-cycling). The global 
dynamics influences the sub-domain evolution by means 
of time- and spatially varying (interpolated) boundary 
conditions. For the details of our AMR algorithm see 
IBarta et all (|2Q10af l. As an illustration, in Fig. [2] the 
regions of enhanced grid resolution at t = 367. Ota (for 
units used see below) are depicted at the background of 
current density and magnetic field. 

The partial differential equations (pQ) are of a mixed 
hyperbolic-parabolic type. We utilized the time-splitting 
approach for their solution: First, hyperbolic (conser- 
vative) part is solved using a second-order FVM leap- 
frog scheme. In a second step the magnetic-diffusivity 
term is solved by means of a (semi -implicit) Al ternating- 
Direction-Implicit (ADI) scheme (|Chungll2002[ ). 

We solve the MHD system of Eqs. (pp) in a 2D simu- 
lation box initially on a global (coarse) Cartesian grid. 
The horizontal and vertical dimensions of calculated box 
are 800 and 6400 grid cells, respectively. Using mirror- 
ing boundary condition at x = in the symmetric CS 
(see below) we obtain a doubled box with an effec tive 
grid of 1600 x 6400 cells (see also IBarta eFaIll2008aL for 
details). We use the following reference frame: The z- 
axis corresponds to the vertical direction, the ?/-axis is 
the invariant (i.e. d/dy = 0) direction along the PIL. 
The x-axis is perpendicular to the current layer and cen- 



tered at the initial current maximum. The simulation is 
thus performed in the xz-plane, while the x?/-plane cor- 
responds to the solar phot osphere; the PIL is located at 
x = 0, z = (see Fig. 5 in IBarta et aTll2008at ). 

The simulation is performed in dimensionless variables. 
They are obtained by the following normalization: The 
spatial coordinates x, and z are expressed in units of 
the current sheet half- width La at the photospheric level 
[z — 0). Time is normalized to the Alfven-wave transit 
time ta = £a/Va,o through the current sheet, where 
Va,o = Bo/^/fiopo is the asymptotic value at x —> oo 
and z = of the Alfven speed at t = 0. Eq. (j2]) for 
anomalous resistivity in the dimension-less variables then 
reads r\ = C(\j\/p — v cr ) for \j\/p > v cr . C and Vqj. are 
now dimension- less parameters. We used C = 0.003 and 
v cr = 15.0 in our simulation. The choice of the threshold 
v cr is not arbitrary as it is closely related to the nu- 
meric resolution reached by the code - see the discussion 
in Section 2J The parameter C was adjusted to reach 
peak anomalous resistivities in the order of 10 6 times 
higher than the Spitzer resistivity in the solar corona - 
similar values for resistivity based on non-linear wave- 
article interaction are indicated by Vlasov simulations 
Biichner fc Elkinall2006l ). 

If not specified otherwise, all quantities in the paper 
are expressed in this dimension-less system of units. In 
order to apply our results to actual solar flares appro- 
priate scaling of dimension- less variables, however, has 
to be performed. The gravity stratification included in 
our model introduces a natural length scale. Assuming 
an ambient coronal temperature of T = 2 MK the corre- 
sponding scale-height for a fully-ionized hydrogen plasma 
is Lq = 120 Mm. The value used in our simulation is 
Lq = 200La, hence La = 600 km. For this scaling the 
flare arcade loop-top is « 10000 km high, which corre- 
sponds well to observed values. The initial CS width 
2La = 1200 km is roughly in line with the fact that the 
CS was formed by stretching of the magnetic field in the 
trail of ejected flux-rope/filament w hich itself has typica l 
transversal dimensions ~ 5000 km (|Vrsnak et al.l 12009). 
It also corresponds (by order of magnitude) to estima- 
tions made from observations of thin layers trailing be- 
hind CMEs that a re sometimes in t erpreted as sign atures 
of current sheets (|Ko et al.ll2003l : iLin et alJl2QQ7h . For 
the ambient magnetic field in the vicinity of the current 
la yer we assume B z q = 40 Gauss (see, e.g., the discussion 
in lKliem et aTll2000l ). 

The initial state has been chosen in the form of a ver- 
tical generalized Harris-type CS with the magnetic field 
B = V X A -\- e qi B qi slightly decreasing with height z 
(IBarta et alJl2Q10ah : 

/ X X \ 

A(x, y, z; t = 0) = -B z0 In f e w cs( z ) + e w cs( z ) j e y 

B y (x,y,z]t = 0) = £3,0(3) 

2; 

p(x,y,z]t = 0) = p exp(- — ) . 

Lg 

In the following we will refer to B x and B z as the prin- 
cipal components and B y as the guide field. The char- 
acteristic width of the initial current sheet varies with z 
as 

d ■ z 2 + z + zq 
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Fig. 4. — Fragmentation of the flare current layer at t = 316. Increasing zoom reveals further smaller magnetic structures (is- 
lands/plasmoids). Red and green lines are positions of B x = 0, and B z = 0, respectively; their intersections represent the X- and 
O-type "null" points. 



and B y o, B z q, po, d, and zq are appropriately chosen con- 



stants: By 



0.2, p = 1.0, B = yjBfo + Bi = 1.0, 

d = 0.003, and zq = 20.0. The initial state given by 
Eq. ((3]) corresponds to a stratified atmosphere in the 
presence of gravity (which is consistent with Eqs. JTJ). 
The divergence of the magnetic field-lines towards the 
upper corona is in agreement with the expansion of the 
coronal field. It also favors up-ward motion of sec- 
ondary plasmoids f ormed in the course of CS tearing 
(Bar ta et al.ll2008bl ). This leads to further filamentation 
of the current sheets which dev elop between the plas- 
moids ([Shibata fc Tanumal l200lh . The current density 
and magnetic field at the initial state are displayed in 
Fig. [3l A rather thick current layer is visible. An en- 
hanced view of the selected area is presented in the right 
panel for the sake of direct comparison with the current- 
density filamentation developed in later stages of the evo- 
lution (see below in Fig. |4]). 

Free boundary conditions are applied to the upper and 
right part of the actually calculated right half of the box. 
It means that von Neumann prescription d/dn = has 
to be fulfilled for all calculated quantities except the nor- 
mal component of magnetic field B n and its contribution 
to the total energy density U. B n and U are extended in 
the second step fulfilling V • B = 0. In order to satisfy 
MHD boundary conditions symmetric (q(—z) = q(z)) 
relations are used for p B y , B Z: U, and anti-symmetric 
(q(—z) = —q(z)) for B x at the bottom boundary while 
velocities are set to zero there (u = 0). This ensures 
that the principal magnetic-field component is vertical at 
the bottom boundary and that the total magnetic flux 
passing through that boundary does not change on the 
rather short time-scales of the eruption, as enforced by 
the pr esence of a dense solar photosphere ([Barta et al.l 
2008a). Mirroring boundary conditions (symmetric in 
p, u y B x , B yi and U and antisymmetric in u x and B z ) 
are used for the left part of boundary at x = (=the 
center of the CS). We use these symmetries to construct 
(mirror) the left half of the full effective box. 

The asymptotic plasma beta parameter at x — >• oo and 
z = 0is/3 = 0.1 and the ratio of specific heats is 7 = 5/3 



(adiabatic response). 

The coarse- mesh sizes are Ax = Az = 0.045 in the 
dimensionless units. Thus, with the reference frame es- 
tablished above the entire box corresponds to (—36, 36) x 
(0, 288) in the xz-plane. The simulation was performed 
over 400 normalized Alfven times. To save disk space 
only the most interesting interval t = 300 — 400 has 
been recorded with a step of 0.5 ta- 

Note that the initial state described by Eq. (j3j) is not 
an exact MHD equilibrium. Nevertheless, the resulting 
field variations are much weaker than those introduced 
by reconnection. 

At the very beginning, in order to trigger reconnection, 
the system is perturbed by enhanced resistivity localized 
in a small region surrounding a line x = 0, z = 30 in the 
invariant direction y for a short time < t < 10 (see also 
iMagara et al.lll996h . This short perturbation sets a local- 
ized inflow which somewhat compresses the current layer 
around the selected point. It should mimic the effect of 
various irregularities that can be expected du ring; the CS 
stretching in actual solar eruptions, see also iRilev et~aTl 
(2007J). Later, the resistivity is switched on only if the 
threshold according to Equation (j2j) is exceeded. As the 
threshold for anomalous resistivity cannot be reached 
for a coarse grid (the threshold for mesh-refinement is 
reached earlier than the threshold for anomalous resis- 
tivity onset), the condition in Eq. ((2]) is actually checked 
only at the smallest resolved scale, for the large-scale 
dynamics we take r\ — 0. 

3. ANALYSIS OF MODEL RESULTS 

We used the above described numerical code in order to 
study, which mechanisms are involved in the transfer of 
free magnetic energy from large to small scales. Thanks 
to the adaptive mesh we were able to cover scales from 
4.5 x 10 -3 La to w 300 La (the larger size of the simu- 
lation box), i.e. over almost five orders of magnitude. 

The early system evolution can be briefly described 
as follows: After the localized initial resistivity pulse 
a flow pattern sets-up that leads to CS stretching (in 
the z-) and compression (in the ^-direction). Eventu- 
ally, the condition in Eq. (j2j) for anomalous resistivity 
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is reached at the smallest resolved scale and first tear- 
ings occur. Dynamics of the plasmoids formed by the 
tearing process leads to further stretching of CSs inter- 
leaving the mutually separating plasmoids. This leads to 
further generation of tearing. Later, after t ~ 300 the 
smallest magnetic structures yet consistent with the res- 
olution start appearing. Here we present an analysis of 
this more developed stage of cascade. Results are shown 
in Figs. H] and [5] Fig. |4] shows the state of magnetic 
field and current density at t = 316. For better orien- 
tation auxiliary lines are added indicating the locations 
where B x = and B z = 0. Their intersections show 
the positions of O-type and X-type "nulls" - the points, 
where only the guide field remains finite. Areas indi- 
cated by blue line are consecutively zoomed (from left to 
right panels). The left-most panel shows the entire sim- 
ulation box and the right-most corresponds to a zoom 
at the limit of the AMR-refined resolution. The figure 
shows, how parts of the current layer are stressed and 
thinned between separating magnetic islands/plasmoids 
formed by tearing instabilities. The current-layer fila- 
mentation is the most clearly pronounced if one compares 
the zoomed views of the same selected area at the initial 
state (Fig.[3j right panel) and the system state at t = 316 
(Fig. Ill the third panel). During the dynamic evolution 
the even thinner, stretched current layers become, after 
some time, unstable to the next level of tearing and even 
smaller plasmoids are formed. The zoomed figures show 
that cascading reconnection has formed plasmoids at the 
smallest resolved scales: The x-sizes of the largest and 
smallest resolved plasmoid in Fig. |4] range from w 10 
down to « 0.01, the z-sizes are from w 0.2 to « 70. 

Plasmoids formed by tearing instability are not only 
subjected to the separation but they can also approach 
each other. As a result the magnetic flux piles- up and 
transversal (i.e. horizontal, perpendicular to the orig- 
inal current layer) current sheets are formed between 



pairs of plasmoids approaching each other. Earlier sim- 
ulations with lower effective resolution treated the plas- 
moid merging as a coalescence process without any in- 
ternal structure of the small-scale (sub-grid) current 
sheet between the magnetic islands s i nce their thickness 
was not resolved dTajima et al .1119871 : iKliem et al.ll200Ci 
iBarta et al.l l20Q8bh . If resolved, however, the transver- 
sal current sheet does not just dissipate. Instead it is 
subjected to the tearing instability in the direction per- 
pendicular to the primary current layer. This is shown in 
Fig. [5l The most detailed resolution (right-most panel) 
clearly reveals the formation of the O-point at x = OLa, 
z = 12.79La and two adjacent X-points at x = 0.37La, 
z = 12.78L A and x = -0.37L A , z = 12.78L A . We 
call this process "fragmenting coalescence" in order to 
emphasize that even smaller structures are formed dur- 
ing the merging of two plasmoids. Thus both the tear- 
ing and (fragmenting) coalescence processes contribute 
to the fragmentation of the original thick and smooth 
current layer. 

3.1. Fragmentation of CS: Scaling 

In order to study scaling properties of the continued 
fragmentation of the magnetic structures associated with 
the current layer we performed both a ID Fourier and a 
wavelet analysis of the magnetic field along the verti- 
cal axis {[x = 0,y = 0,z G (0,288)]}. We use for this 
study the B x component since B z = there due to the 
boundary condition. The results are shown in Fig. [6l 
The upper panel shows the magnetic field and current 
density in the sub-set of the entire computation domain 
(z G (0,200), note the rotated view), where the current 
layer is fragmented. Panel (b) shows the profile of B x 
along the current-layer axis, and panels (c) and (d) the 
Fourier and wavelet analyses of this profile. The Fourier 
power spectrum exhibits a power-law scaling with the 
spectral index s = —2.14 in rather wide range of scales 
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Fig. 6. — The ID scale analysis of the magnetic field structure 
along the line x = 0. (a) Magnetic field lines and the current den- 
sity structure (green) in the computational domain at t = 316. The 
z-axis shows positions both in the units of La (top) and in kilome- 
ters according to scaling adopted in Section [5| (b) Profile of the 
B x component of magnetic field along the line x = 0. (c) Fourier 
power spectrum of the B x profile, (d) Wavelet power spectrum of 
the B x profile. 



300 km - 10000 km. This clearly indicates cascading 
nature of the continued fragmentation. 

The energy-transfer cascade ends at w 300 km in 
Fig. [6] This is closely related to the dissipation threshold 
v cr that has been chosen as v cr = 15.0 in our simulation. 
By selecting this value we shifted dissipation-scale do- 
main into the window of resolved scales - see the discus- 
sion in Section HI Typical width of dissipative current 
sheets in our model is thus « L^/15.0 = 40 km. 
Since the plasmoid dimensions along the CS are about 
one order of magnitude (typically 6 x ) larger than across 
CS, distribution of magnetic energy in structures along 
CS, which is depicted in Fig. [6j reaches its dissipation 
scale at ~ 6 x 40 km = 240 km. In reality, the ion iner- 
tial length d{ = cjuj vi is co nsidered as a t ypical width of 
dissipative current sheets (Buchner 2007). Its value for 




x/L^ x/L^ 



Fig. 7. — Filament at ion and splitting of diffusive regions at 
t = 328. Detailed treatment shows how the magnetic diffusivity is 
concentrated into multiple thin channels. The blue-line bounded 
areas indicate the diffusion regions where the generalized anoma- 
lous resistivity (Eq. [2} is finite. Field lines are dashed to be better 
distinguished from the resistive region boundaries. Current den- 
sity is represented by gray levels. The right panel displays detailed 
view of the rectangular box in left. Red and green lines are as in 
Fig. HI 



parameters Bo and Va used in this paper is di ~ 4 m in 
the simulation-box center, i.e. at z ~ 140 (see Fig. [I]). 

The most pronounced features in the wavelet spectrum 
are the locations of low signal (the white islands). They 
correspond to the filamented parts of the current layer 
between plasmoids. Their distribution indicates, that the 
filamented current sheets are embedded within the global 
current layer in a hierarchical (qualitatively self-similar) 
manner. 

3.2. Fragmentation of CS: Diffusion regions 

The current sheets are filamented down to the reso- 
lution limit of our simulations (in reality to the kinetic 
scales). The smallest current-density structures contain 
dissipative/ acceleration regions. In the following we will 
study the structure, distribution and dynamics of these 
non-ideal regions embedded in the global current layer. 

Cascading reconnection and consequent fragmentation 
of the current layer may have significant impact also for 
particle acceleration in solar flares. Instead of a single 
diffusion region assumed in the 'classical' picture of the 
solar reconnection, cascading fragmentation causes the 
formation of large amount of thin non-ideal channels. 
The structuring of non-ideal regions in our simulation is 
depicted in Fig. [71 The left panel shows two areas of 
dissipation around x = 0, z = 36 and x = 0, z = 44. A 
closer look (right panel, note the large zoom), however, 
reveals that the bottom dissipation region is structured 
and it is in fact formed by two regions of finite mag- 
netic diffusivity that are associated with two X-points 
at x = 0, z = 35.40 and x = 0, z = 35.97 interleaved 
with a (micro) plasmoid. The multiple dissipative re- 
gions embedded in the global current layer are favorable 
for efficient (and possibly multi-step) particle accelera- 
tion. At the same time they provide a natural expla- 
nation of fragmented energy release as it has been in- 
ferred from HXR and rad io observations ( Asch wandenl 
120021 : iKarlickv et al.|[2000f ). Since they are embedded in 
the large-scale current layer the 'classical' well-organized 
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global picture of eruptions is kept simultaneously. 

Fig. [7] shows that the X-points formed in the thinned 
current sheets between magnetic islands are connected 
with the thin channels of magnetic diffusivity. Hence it 
is appropriate to study the distribution and dynamics 
of these non-ideal regions by means of tracking the X- 
points associated with them. We present such analysis in 
Fig. [51 In order to see the "skeleton" of the reconnection 
dynamics we followed the positions of all magnetic "null" 
points during the entire recorded interval t = 300 — 400. 
The results show the kinematics of the O-type (red cir- 
cles) and X-type points. Motivated by our endeavour 
to establish a closer relation of the model to observable 
quantities in our consecutive study ([Barta et al.ll2010bl ) 
we also paid special attention to the magnetic connec- 
tivity of the X-points to the bottom boundary. For this 
sake, the X-points connected to the model base (= the 
photosphere) are painted as green asterisks while the un- 
connected X-points are displayed as black crosses. 

Since we are interested in the cascade that already 
developed into the stage reaching the smallest resolved 
structures (t > 300), many X- and O- points connected 
with the larger-scale (and therefore longer-living) struc- 
tures (plasmoids) are already formed and their space- 
time trajectories enter Fig. |H(a) from the bottom. Fig. [8] 
thus shows mainly full life-cycles of the X- and O- points 
related to the smallest resolved plasmoids. It is best 
visible in the three bottom panels (b) - (d) that show 
zoomed views (projected to the zt-plane) of typical exam- 
ples of null-point dynamics - the creation of temporary 
X-0 pairs (panels (b) and (d)) and plasmoid merging (c). 
As it can be seen from panels (b) and (d) the X-points 
can become magnetically connected (the right X-point in 
panel (d)) or disconnected (panel (b)) to/from the bot- 
tom boundary during their lifetime. Note also the split- 
ting (and subsequent merging) of the X-point at x = 0, 
z = 12.8 into an X-O-X configuration between t « 360 
and t ~ 380 in panel (a). This process maps the tearing 
in the transversal (horizontal) current sheet formed be- 
tween interacting plasmoid and the loop-arcade (see also 
Fi g;. [5]). Note that FigJS] can be compared with Fig. 5 
in ISamtanev et all ((2009). The main difference is just 
in the presence of the off-plane X-points formed by the 
fragmentation of the CS between coalescing plasmoids in 
our simulation. 

4. DISCUSSION 

Reconnection in the trailing current layer behind an 
ejected flux-rope (filament) is a key feature of the 'stan- 
dard' CSHKP scenario of solar flares. A large amount of 
free magnetic energy is accumulated around this rather 
thick (relative to plasma kinetic scales) and very long 
layer. The thickness of this l ayer was estimated both 
from the observed brightening (|Ko et al.ll2003l : iLin et aD 
120071 ) an d based on a typica l transversal dimensions of a 
filament (|Vrsnak et al.ll2009l ). Both ways one obtains the 
order of magnitude of ~ 1000 km. On the other hand col- 
lisionless reconnection requires dissipation at very small 
scale, thin current sheets with typical width o f the order 
of ~ 10 m in the solar corona (Buchner 2007). The fun- 
damental question arises how the accumulated energy is 
transferred from large to small scales. Or, in other words, 
what are the mechanisms of direct energy cascade in 
magnetic reconnection. We addressed this question using 



high-resolution AMR simulation covering broad range of 
scales to investigate the MHD dynamics of an expanding 
current layer in the solar corona. 

4.1. Mechanisms of direct energy cascade 

Our simulations reveal the importance of a contin- 
ued fragmentation of the current layer due to the in- 
teraction of two basic processes: The tearing instabil- 
ity of stretched current sheets and the fragmenting co- 
alescence of flux-ropes/plasmoids formed by the tear- 
ing and subsequently forced to merge by the tension of 
ambient magnetic field. After ejection of the primary 
flux-rope (i.e. the filament/CME), a trailing current 
layer is formed behind it which becomes long and thins 
down. As it has been pointe d out by theoretical anal- 
ysis by lLoureiro et al.l (|2007[ ) , current layers with high 
enough length-to-width ratio become unstable for fast 
plasmoid instability. Moreover, any irregularity in the 
plas ma inflow that stretches the sheet facilitates the tear- 
ing ([Lazarian fc Vishniacl [l999). 

Plasmoids that are formed are subjected to the ten- 
sion o f ambient magneti c field, which causes them to 
move ([Barta et al.ll200 8b). The motion can lead to their 
increasing separation. A secondary current layer then 
formed between them becomes, again, stretched and a 
secondary tearing instability can take place. This simu- 
lation result, illust rated by Fig.lH fully confir ms the sce- 
nario sugge sted by Shibata & Tanuma (2Q0 l|), developed 
furthe r by lLoureiro et al.l (|2007l ) and lUzdenskv et al.l 
(2010| into the analytical theory of chain plasmoid in- 
stability. The results are also in qualitative agree- 
ment with the simula tions of plasmoid instab il ity b y 
ISamtanev et al.l (|2QQ9h and iBhattacharjee et all (|2QQ9h . 
which has been performed, however, with constant resis- 
tivity. 

In addition to that, our simulation has shown that the 
converging motion of plasmoids leads to a magnetic-flux 
pile- up between mutually approaching plasmoids. Con- 
sequently, secondary (oppositely directed) current sheets 
are formed perpendicular to the original current layer. 
While previous studies found only unstructured current 
density pile-ups between merging magnetic islands our 
enhanced-by-AMR resolution reveals secondary tearing 
mode instabilities that take place in the transversal to the 
primary current sheet direction (see Fig. [5]). This pro- 
cess represents a new mechanism of fragmentation and 
changes our view to coalescence instability, which has 
been hitherto commonly considered as a simple merg- 
ing process of two plasmoids contributing to the inverse 
energy cascade only. Note that this behavior is differ- 
ent from that seen f or plasmoids at the dissipation scales 
in PI C simulations (Drake et al. 2005; Karlick v fc Barta] 
120071 ). where plasmoids merge without subsequent tear- 
ing. 

One can suppose that with even higher spatial resolu- 
tion one would see more subsequent tearing mode insta- 
bilities altering with the fragmenting coalescence of the 
resulting magnetic islands/flux-ropes. As a result third- 
and higher order current sheets could form. 

To sum up, the results of our simulation support 
the i dea that both the t earing ([S hibata fc Tanuma] 
120011 : lLoureiro et "aD 120071 : lUzdenskv et al.l I2010D and 
"fragmenting coalescence" processes lead to the for- 
mation of consecutively smaller magnetic structures 



9 




z/L A z/L A z/L A 



Fig. 8. — Kinematics of the magnetic "null" points (X-points/dissipative regions and O-points/plasmoids). Red circles denote O-point 
positions, black crosses the X-point positions and green asterisks show the positions of those X-points which are magnetically connected 
to the photosphere. Three bottom panels show detailed view of selected rectangles and represent typical processes of X- and O-point 
dynamics: Bifurcation (null-pair creation) and merging (null-pair annihilation). 




Incompressible FD: Sheared Flow 
rotv = £2 divv=0 



MHD: Current Layer 
rotB = |i i divB = 



Fig. 9. — Schematic view of cascade of fragmentation of the cur- 
rent layer by tearing and driven coalescence processes. Increasing 
zoom shows similar kind of processes repeating on smaller scales. 
Red and blue areas of various hue indicate various intensities of 
positive and negative out-of-plane (j y ) component of current den- 
sity. 



(plasmoids/flux-ropes) and associated current filaments. 
Subsequent stretching and compression cause a filamen- 
tation of the current. This situation is schematically de- 



Cascade of vortices. 



Cascade of plasmoids? 



Fig. 10. — Illustration of analogy between sheared-flow in in- 
compressible fluid dynamics (FD) and magnetic reconnection in a 
large-scale current layers. In the case of FD the mechanism of en- 
ergy transfer from macroscopic to dissipative (molecular) scale is 
known - it is the cascade of vortex-tubes. The cascade of magnetic 
flux-ropes/plasmoids can play this role in the case of magnetic dis- 
sipation. 



picted in Fig. [9] which can be seen as a genera l izatio n 
of the scheme in Fig. 6. in Shib ata fc Tanumal (|200l[ ). 
One can expect that this cascade will continue down to 
the scales where the magnetic energy is, finally, dissi- 
pated. Note that the physics and the corresponding scal- 
ing laws may change at intermediate (but still relatively 
small) scales when additional contributions to the gener- 
alized Ohm's law becom e significant, e.g. a Hall te rm - 
see recent simulatio ns bv lShepherd fc Cass ak (2010) and 
iHuang et all (|2010l ). 

4.2. Impact on reconnection efficiency 



10 



Reconnection in the current sheet between merging 
plasmoids is fast since it is driven by ambient-field mag- 
netic tension which naturally pushes the flux-ropes to- 
gether. Thus even shorter time-scales can be reached 
by t his process than by tearin g cascade in the stretched 
CS (|Shibata fc Tanumal l200lh . And yet another point 
makes the overall reconnection process more efficient: 
Many magnetic-flux elements - except those ejected out- 
wards to the escaping CME - reconnect several times. 
First, during the primary tearing and plasmoid forma- 
tion and then again during plasmoid coalescence. Since 
coalescence leads to a follow-up tearing instability (in 
the transversal direction) the remaining magnetic flux is 
subjected to another act of magnetic reconnection. This 
process rese mbles the recurrent separator reconnection 
simulated by iParnell et al.l (|2QQ8f ) . 

4.3. Relation to turbulence onset 

To some extent the initial situation of global, smooth 
and relatively thick sheets is similar to the turbulence 
on-set in a sheared flow in (incompressible) fluid dynam- 
ics (FD) as schematically shown in Fig. [TUl Usually the 
typical length-scale of shear flows - the counterpart of 
the width of current layers - is much larger than the 
dissipative (molecular) scale. The mechanism of energy 
transfer from large to small scales in classical FD is me- 
diated by a cascade of vortex tubes: Large-scale vor- 
tices formed by shear flows can mutually interact giving 
rise to increased velocity shear at the smaller scales in 
the space between them. Each small shear flow element 
formed by this process can be, again, subjected to this 
fragmentation. Based on our simulation results, we sug- 
gest a similar scenario for current-layer fragmentation. 
The role of the vortex-tubes in FD is in MHD taken over 
by flux-ropes/plasmoids. 

In analogy with the on-set of turbulence in sheared 
flows, one could expect that a dynamical balance would 
arise between fragmentation and coalescence processes in 
later more developed stage. This should be manifested 
by a power-law scaling rule. Using AMR we reached a 
rather broad (five orders of magnitude) range of scales. 
This allowed us to perform a ID scaling analysis of the 
magnetic-field structures formed along the current layer 
for the first time. The scaling rule found exhibits, in- 
deed, a power- law distribution with the index s = —2.14 
(Fig. [6}. Since our resolution still does not allow to make 
this scaling analysis only within the small selected sub- 
domain around the CS center where one could expect 
isotropic 'turbulence' (we would lack sufficient range of 
scales for that) it is difficult to compare the spectral in- 
dex found over the whole (clearly anisotropic) simulation 
domain with the values expected from the theory for fully 
developed isotropic turbulence. 

Obtained power-law distribution is also in qualitative 
agreement with the conc ept of fractal reconnection by 
Shib ata fc Tanumal ([2001) and wit h hierarchical analyt- 
ical model of pl asmoid instabili t y dLqureiro et aLll2007l ) 
as described by lUzdenskv et al.l (|2010l ). Thev use distri- 
bution functions for plasmoid width and contained flux 
in order to characterize statistical properties of plasmoid 
hierarchy rather than power spectrum. We plan to per- 
form similar analysis of our results in the future study in 
order to compare the results also quantitatively. 

In order to obtain as broad as possible scale range 



in the plane where reconnection occurs, we performed 
these simulations using 2.5D approach. The question 
arises to what extent a full-3D treatment would change 
the resulting picture. In the FD cascade vortices are 
deformed, their cross-sections change along their main 
axis, even in the topological sense. The object de- 
fined as a single vortex tube in one place can be split 
into two in another location. One can expect a sim- 
ilar behavior of plasmoids/flux-tubes in MHD. There 
they could be subjected to the kink and similar in- 
stabilities with k y > 0. Such processes would natu- 
rally lead to the modulation of reconnection rate along 
the PIL. Observations indicating such effect ha ve al- 
ready been presented (|McKenzie fc Savagd l20Q9h . To 
some extent the expected behavior can also be ob- 
tained for kink instabilities of tiny current channels at 
the dissipative scale studied by 3D PIC sim ulations 
(IZhu fc Windeel 119961: iKarlickv fc Bartal l2008h. Never- 
theless. lEdmondsonet al.l (|2010l ) found no such evidence 
in their 3D AMR MHD simulation, which uses, however, 
a different set-up. Sizes of the plasmoids formed under 
3D perturbation in the direction that corresponds to the 
invariant ?/-axis in our 2.5D case have been found very 
short preventing a kink-like instabilities to develop. On 
the other hand, the resulting plas moid lengths might de- 
pend on the initial g uide field ([Edmondson et all 120101 : 
iDahlburg etllll2005l i. 

To sum up, the answer to this question can be found 
only via full 3D simulations with similar initial set-up as 
we used here. Therefore we plan to extend our current 
2.5D simulations with very high in-plane resolution with 
moderately resolved structuring in the third dimension. 

4.4. Fragmented energy release and particle acceleration 

Cascading fragmentation of the current layer is closely 
related to another puzzling question of current solar 
flare research - the apparent contradiction between ob- 
served regular large-scale dynamics and signatures of 
fragmented energy release in (eruptive) flares. This du- 
ality is reflected by two classes of flare models: The 'clas- 
sical' CSHKP scenario based on magnetic reconnection 
in a single global flare current sheet and the class of 
"self-organized criticality" (SOC) models based on the 
avalanche of small-scale reconnection events in multiple 
current sheets formed as a cons equence of either chaotic 
(|Aschwan den 20021: IVlafiosll2007l ) or regular but still com- 
plex boundary motio ns causing, e.g., magnetic braiding 
(|Wilmot-Smith et al.ll20Toh . 

The model of cascading reconnection has the poten- 
tial to provide a unified view on these seemingly very 
different (see the discussion in the next paragraph) ap- 
proaches. From the global point of view, it coincides 
with the classical CSHKP model keeping the regular pic- 
ture of the process at large scales. At the same time, 
due to the internal current-layer fragmentation the tear- 
ing/coalescence cascade forms multiple small-scale cur- 
rent sheets and potential diffusive regions. As a conse- 
quence, fragmented energy release, e.g., by particle accel- 
eration, can take place in these tiny regions. To some ex- 
tent this finding can be seen a s a follow-up of the 'bursty' 
reconnection regime found bv lKliem et al.l (|200Q[ ). These 
authors show that intermittent signal (X-ray, radio) can 
be related to chaotic pulses of the (resistive) electric field 
in the dissipation region around X-point. The pulsed 
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regime is a consequence of non-linear interplay between 
governing MHD equations and the anomalous resistiv- 
ity model. In our view, however, these pulses result from 
the subgrid physics unresolved in earlier simulations. Es- 
sentially, what has been seen as a single dissipative re- 
gion around a single X-point in the coarse-grid models 
is in fact a ("fractal-like") set of non-ideal areas around 
multiple X-points (see Fig. [7j) that interleave very small- 
scale mutually interacting plasmoids. This view based on 
high-resolved simulation is perhaps closer to the term of 
"fragmented energy release" that assumes the energy dis- 
sipation to be performed via many concurrent small-scale 
events appearing in multiple sites distributed in space. 
In this context it is interesting to note how surprisingly 
we ll the phenomenolo gical resistivity model (Eq. [2j) used 
bv lKliem et al.l ((2000) mimics the sub-grid scale physics 
as it has been able to reproduce qualitatively temporal 
behavior of resistive electric field even without resolving 
actual processes that are responsible for it. Note also 
that a possible role of tearing and coalescence in frag- 
mentation of the ener gy release in solar flares has been 
mentioned already by iKlieml ([ 19901 ) . 

We would like to emphasize that there is a fundamen- 
tal difference between the fragmented energy release by 
cascading reconnection and SOC models. It is rooted 
in the fact that in cascading reconnection the complex- 
ity /chaoticity is due to an intrinsic current-layer dynam- 
ics, i.e., due to spontaneous fragmentation, while it is 
introduced in SOC models through (external) boundary 
conditions (chaotic boundary motions). In fact, these 
two concepts are contrary to some extent: While in 
SOC-based models the global flare picture is built as 
an avalanche of many small-scale events (bottom— >-top 
process in the scale hierarchy) in cascading reconnection 
small scale structures are formed as a consequence of in- 
ternal dynamics of large-scale CS (top^bottom process). 

Fragmented energy release is closely related to the 
number problem of particles accelerated in solar flares. 
A single diffusion region assumed in the CSHKP model 
provides a far too small volume for accelerating strong 
fluxes of particles as they are inferred from the HXR ob- 
servations. This argument has been used in favor of SOC- 
based models as they provide energetic-particle spectra 
and time-profile distributions as observed and explain 
large energetic particle fluxes. 

We suggest that, however, the inclusion of cascading 
reconnection into the CSHKP has even more capabilities 
than SOC models. It could explain both the distribu- 
tion and the number of accelerated particles, based on a 
physical consideration of many small-scale current sheets 
which can host tiny diffusive channels that all can act as 
the acceleration regions (see Fig. [7J). 

Here it is appropriate to make one technical comment: 
In an MHD simulation with resistivity model described 
by Eq. (J2j), respectively its dimension- less version, the 
size and the number of diffusive regions are controlled 
mostly by the threshold v cr for the onset of (anomalous) 
diffusivity. The higher v cr is chosen, the thinner the cur- 
rent sheets can become, the smaller and more numerous 
are the embedded diffusion regions. Since one has to re- 
solve these diffusive regions in the simulation, one has to 



choose the threshold v cr low enough to be able to resolve 
the dissipation regions appropriately. In ideally resolved 
simulations, covering all scales down to the real physi- 
cal dissipation length, the critical velocity v cr could be 
chosen of the order of phy sically releva nt value - the elec- 
tron thermal speed vr e (|Buchnerll2QQ7l ). In dimensionless 
units this corresponds to v cr = L\/di \fmi/m e \/ 'ft /2, 
where ra$ and m e are the proton and electron masses. For 
a technically limited spatial resolution, one has to choose 
a (much) smaller value of v cr in order to resolve the small- 
est possible current sheets, before dissipation sets in, by 
a reasonable number of grid points. Since the resolution 
in our current AMR simulation is higher than in ear- 
lier models, we could choose a more reasonable value of 
v cr = 15.0 (while older simulations used y cr — 3.0 - see, 
e.g. JKliem et aLll2QQ0t fearta et al.ll2008at ). This allowed 
us to track down more fragmented, smaller reconnection 
regions. If we extrapolate this trend, we can expect that 
with even higher resolution one would find even more 
and tinier diffusive regions. They would be grouped hi- 
erarchically (self-similar ly), occupying a sub-space of the 
global current layer. Such kind of distribution is indi- 
cated in the wavelet spectra (white islands in Fig. [6fd)), 
and also by the positions and motion of the associated X- 
points in Fig. [51 The latter shows a structured grouping 
of "null points" and their various life times. 

5. CONCLUSIONS 

Our simulation has shown that cascading reconnec- 
tion due to the formation and fragmenting coalescence of 
plasmoids/flux-ropes is a viable physical model of frag- 
mented magnetic energy release in large-scale systems, 
like solar flares. Cascading reconnection addresses at 
once three key problems of the current solar-flare re- 
search: The scale-gap between energy-accumulation and 
dissipation scales, the duality between regular global- 
scale dynamics and fragmented energy-release signatures 
observed simultaneously in solar flares, and the issue of 
particle acceleration. All these problems arising from 
observations appear to be tightly related via cascading 
reconnection. 

In order to evaluate relevance of the cascading recon- 
nection for actual solar flares further it is desirable, how- 
ever, to identify and predict model-specific observables 
and to search for them in observed data. We are going 
to propose possible specific signatures and co mpare them 
with o bservations in a consecutive paper ([Barta et al.l 

HoToS). 
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